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摘要 :[ 目 的 ] 实践 中 95% 以 上 的 工业 过 程控 制 问 题 都 可 以 通过 PID 控制 算法 解决 . 本 文 在 继 
承 并 创新 地 运用 二 自由 度 PID 内 模 控 制 技术 的 基础 上 , 建立 了 一 种 通用 控制 器 优化 整定 新 方 
法 . [方法 ] 方 法 中 引入 模范 系统 , 离线 仿真 该 模范 系统 后 获得 优化 模板 ; 在 信号 激励 实际 控制 系 
统 的 过 程 中 , 用 该 模板 引导 特定 算法 整定 PID 内 模 控 制 参数 . [结果 ] 在 无 精准 对 象 参 数 的 情况 
下 , 经 数 次 循环 迭代 后 , 可 同步 地 实现 对 象 参 数 的 辨识 与 控制 性 能 优化 . [结论 ] 此 法 整定 效率 
高 , 便于 后 期 系统 在 线 维护 , 减少 了 对 实施 人 员 技 术 与 经 验 的 要 求 . 所 用 控制 器 继承 内 模 控 制 
大 时 潇 控制 效果 好 、 重 棒 性 强 的 特点 , 兼顾 ”目标 跟踪 ”与 ”干扰 抑制 ”两 种 性 能 优化 . 控制 算 
法 结构 简单 、 直 观 , 易于 在 原 PID 控制 系统 升级 改造 或 用 嵌入 式 系 统 软 、 硬 件 实现 , 便于 在 生 
产 中 推广 应 用 .[ 局 限 与 改进 方向 ] 受 限于 传统 内 模 控 制 的 设计 要 求 ， 对 于 具有 负 响 应 特性 的 受 
控 对 象 , 该 方法 需 采 用 预 补偿 处 理 机 制 方 能 适用 或 者 寻找 新 的 模范 系统 模型 并 研究 改进 算法 . 
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Abstract: [Objective]Imn practice, more than 95% of industrial process control problems can be 
solved by PID control algorithm. In this paper, a new method of general controller optimization 
setting is established on the basis of inheriting and innovatively using 2DOF PID internal model con- 
trol technology. [Methods]Model System was introduced into the method, and optimization template 
was obtained after offline simulation of the Model System. In the process of the signal excitation of 
the actual control system, the PID internal model control parameters are set by a specific algorith- 
m guided by the template. [Results]Without accurate object parameters, the identification of object 
parameters and optimization of control performance can be completed simultaneously after several 
cycles of iteration. [Conclusions]This method has high setting efficiency, is convenient for later on- 
line maintenance of the system, and reduces the requirements for the technology and experience of 
the implementation personnel. The controller inherits the characteristics of good control with large 
time delay and strong robustness of internal model control, and both ”target tracking” and ”inter- 
ference suppression” performance optimization. The control algorithm is simple, intuitive, and easy 
to be upgraded in the original PID control or realized by embedded system software and hardware, 
which is easy to be popularized and applied in production. [Limits and future work]Restricted by 
the design requirements of traditional internal model control, for the controlled object with negative 
response characteristics, this method should adopt the pre compensation processing mechanism or 
find a new model system and study the improved algorithm. 
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1 引言 

牛顿 第 二 定律 在 系统 建 模 中 起 关键 作用 , 其 对 应 的 动力 学 模型 通常 为 一 、 二 阶 微分 方程 . 在 工业 生产 过 程 中 
常见 的 也 多 是 一 、 二 阶 惯性 加 时 滞 环 节 (FOPDT、SOPDT) 受 控 对 象 , 因此 实践 中 95% 以 上 的 过 程控 制 问题 都 可 
以 通过 PID 的 控制 算法 解决 上 1. 人 们 对 PID 控制 器 进行 了 大 量 的 研究 71, 但 是 迄今 为 止 人 们 对 PID 控制 的 基 
础 理论 和 参数 设计 方法 研究 远 未 成 熟 . PID 控制 器 实现 中 最 具 挑战 性 的 任务 之 一 是 如 何 设计 控制 器 三 个 参数 : 
PID 控制 中 比例 项 主要 用 来 镇 定 系统 , 积分 项 可 消除 系统 静态 偏差 , 微分 项 可 预测 未 来 误差 趋势 . 因 其 综合 的 控 
制 能 力 , 使 PID 控制 器 恰好 适用 于 上 述 受 控 对 象 的 系统 控制 , 并 对 局 部 非 线性 对 象 也 适用 . 这 三 个 参数 传统 上 是 
通过 实验 、 经 验 积累 或 者 两 者 结合 进行 调 优 的 , 在 对 象 模型 参数 未 知 的 情况 下 , 多 依赖 整定 人 员 的 经 验 在 现场 
十 几 次 甚至 几 十 次 地 反复 调整 优化 控制 参数 ; 或 者 通过 估 测 有 具体 对 象 的 模型 后 , 用 各 种 经 验 公 式 或 者 模拟 仿真 
方法 对 控制 参数 不 断 优化 选 定 & 10. 其 控制 效果 由 于 实施 人 员 的 技术 、 经 验 或 模型 估 测 的 精度 等 多 种 因素 影响 
而 差别 很 大 , 且 整 定 效 率 不 高 . 

Garcia 和 Meorari 在 1982 年 提出 一 种 控制 器 优化 设计 方法 : 内 模 控 制 11, 是 一 种 实用 性 很 强 的 控制 方法 . 其 结 
构 简 单 、 设 计 直 观 简便 、 调 节 参 数 少 , 特别 是 鲁 棒 性 、 抗 干扰 性 能 以 及 大 时 灌 系 统 控制 效果 的 改善 尤为 显著 . 
但 缺点 是 不 能 兼顾 ” 目标 值 跟踪 ”和 ”扰动 抑制 ”两 种 最 佳 控制 性 能 . I. M. Horowitz 在 1963 年 提出 了 二 自由 度 
结构 PID 的 概念 05. 采用 二 自由 度 PID 结构 的 内 模 控 制 (简称 IMC-2DOF-PID) 进行 控制 器 设计 是 一 大 进 
步 03-19, IMC-2DOF-PID 不 仅 具 有 内 模 控 制 器 的 优点 , 也 兼顾 了 上 述 两 种 性 能 的 优化 ; 并 有 利于 对 原 有 的 PID 
控制 系统 进行 改造 , 也 便于 用 骸 入 式 系统 软 、 硬 件 实现 . 因此 为 广大 工程 人 员 所 接受 05-19. 然而 实现 此 方法 的 
前 提 是 必须 获取 精准 的 对 象 模 型 参数 . 

控制 技术 发 展 到 今天 , 其 理论 与 应 用 取得 了 长 足 的 进步 , 涌现 出 许多 理念 先进 、 手 段 多 样 、 算 法 智能 的 各 种 
新 型 控制 方法 上 9. 只 是 这 些 方法 在 大 范围 推广 到 实际 生产 时 还 是 遇 到 很 大 的 阻碍 , 原因 或 是 需 精准 辨识 对 象 参 
数 模型 , 或 是 要 求 有 丰富 的 实践 经 验 与 专业 技能 的 实施 人 员 在 现场 逐个 针对 每 个 受 控 对 象 反复 测试 与 调整 才能 
得 到 可 用 的 控制 效果 . 由 于 实际 受 控 对 象 在 长 期 运行 中 存在 损耗 , 为 维护 系统 保持 原 有 性 能 仍 需 不 定期 地 重复 
这 些 繁 琐 作 业 . 以 传统 内 模 控 制 的 实现 流程 为 例 : 首先 需要 在 现场 针对 具体 对 象 辨识 其 模型 , 而 后 在 设计 中 根据 
可 实现 要 求 增加 一 些 蔡 代 和 补偿 的 滤波 环节 , 增加 环节 的 参数 还 需要 回 到 现场 结合 实际 对 象 对 系统 进一步 调 优 . 
两 次 现场 实施 过 程 都 需要 有 经 验 的 技术 人 员 才 能 完成 . 人 为 因素 影响 较 大 、 效 率 不 高 , 后 期 维护 也 不 方便 . 

如 能 将 对 象 辨识 和 系统 优化 合并 为 一 个 过 程 在 线 上 系统 的 闭环 控制 中 施行 , 且 通 过 算法 得 出 优化 的 控制 参 
数 , 就 可 以 明显 提高 工作 效率 , 减少 人 为 因素 影响 , 便于 在 线 维护 系统 , 也 有 利于 实际 应 用 的 推广 . 为 此 目的 作者 
在 继承 二 自由 度 PID 内 模 控 制 设计 技术 的 基础 上 , 建立 一 种 通用 控制 器 优化 新 方法 L131; 方法 中 特别 引入 模范 系 
统 概念 , 优化 其 模范 参数 并 制 成 响应 模板 , 该 响应 模板 适用 于 大 多 数 对 象 . 在 现场 信号 激励 控制 系统 过 程 中 , 借 
助 此 响应 模板 , 由 特定 算法 友 代 地 整定 PID 三 个 控制 参数 , 最 终 同步 实现 控制 优化 和 对 象 参 数 辨识 . 辨识 得 到 对 
象 参数 后 还 可 借助 不 同 的 响应 模板 快速 改变 控制 器 各 种 优化 性 能 , 进一步 提高 实施 效率 . 尽管 该 方法 的 实现 机 
理 和 实施 流程 有 很 大 变化 , 但 控制 器 仍 是 PID 的 结构 形式 , 算法 简单 易于 实现 ; 控制 特性 继承 了 内 模 控 制 的 优势 
特点 . 下 面 将 对 该 方法 的 理论 依据 、 技 术 设 计 思 路 和 推导 过 程 进行 阐述 , 且 通 过 模拟 仿真 验证 其 效果 并 分 析 总 
结 其 特点 : 
2 通用 二 自由 度 PID 内 模 控制 系统 的 设计 

在 工业 生产 过 程 中 存在 着 大 量 型 如 (]) 式 的 一 、 二 阶 惯性 加 时 灌 环 节 的 受 控 对 象 : 

Gu(5) = 二 STARSHRIe 

(KK2>0iKI>0iKo>0ir> 二 0) 


可 用 图 1 结构 的 二 自由 度 PID 控 制 系统 进行 控制 . 


一 
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图 1. 含 前 置 滤波 的 二 自由 度 PID 控制 系统 结构 图 


Fig. 1. Two degrees of freedom PID control system containing pre filtering structure 
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其 控制 器 形 如 (2) 式 , 滤波 器 如 (3) 式 , 系统 响应 如 (4) 式 . 


Ge(S)= EK, (1+ as +7T5)/ (TS+1) 


一 (X25+1) 
G7(5) = C+ 


Y(5) = Gyr sR(S)+ 


男 一 种 含 滤波 反馈 的 内 模 控 制 系统 结构 如 图 2 所 示 . 


D(S) 


1 
EE 


图 2. 滤波 反馈 内 模 控 制 结构 图 


Fig. 2. Filter feedback internal model control Structure 


图 中 内 模 部 分 G,,(5) 等 于 受 控 对 象 G,(5), 形 如 (5) 式 . 

Gm(S) 二 G»(S) 一 Gm-(S)Gm+(S) 
式 中 Gm_(5) 是 对 象 的 最 小 相位 部 分 ; Gm (5S) 是 时 淆 环节 的 Taylor 一 阶 近 似 , 如 (6) 式 . 
为 (7) 式 : 


> 于 

) ss Ko2S52+K1S+Ko 
)=e-" NT 一 79 
) 


2 1—7S 
K2S2+K1S+Ko 


依据 内 模 控制 实现 要 求 , 设计 内 模 控 制 系统 中 的 Q.(5),Qy(S), 形 如 (8), (9) 式 : 


Q (5) ee G(s) _ K2S?+K1S+Ko 
有 (A15+1)(A25+1) (和 1S 二 1)(X2S 十 1) 


二 (X19 二 1) 
9r(9) = sy 


该 内 模 控 制 系统 可 以 等 效 为 图 3 的 结构 . 


图 3. 等 效 为 前 置 滤波 二 自由 度 PID 控制 系统 结构 图 


Fig. 3. Equivalent to pre filter two degrees of freedom PID control system structure 
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(2) 
(3) 
(4) 


(5) 


因此 Gm(S) 近似 


(6) 


(7) 


(8) 
(9) 


图 中 虚线 部 分 整体 可 看 作 一 个 控制 器 Ge( 5), 因此 内 模 控制 系统 和 图 1 的 二 自由 度 PID 控制 系统 的 结构 之 间 就 
有 (10), (11) 式 的 对 应 关系 : 


Gy == 太 (10) 
一 Qf Qe 
Ge = Odom (11) 


将 内 模 G,(5) 及 设计 的 内 模 控制 器 Q.(5), 滤波 器 Qy(S) 〈 即 (7), (8) 和 (9) 式 ) 代入 (11) 式 可 得 到 图 1 控制 系 
统 的 PID 控制 器 如 下 : 


Ge(5) K2S2+K1S+Ko 


S (X35S+2X2+7) 
— _K (1+RS+RS) 
32H7F ~ (12) 
呈 下 元 
(二 s+t1) 


对 照 (2) 式 和 (12) 式 的 形式 可 得 PID a 
Ke.= zs Zi (K1 > 0) 
TT = Ko/Ki;(Ko > 0;K >0) 
Ty = Kz/Ki; (Ki > 0; K, > 0) 
入 
Tf = 双生 
当 受 控 对 象 为 积分 型 (K0 = 0) 时 , 该 控制 器 无 积分 控制 项 ; 如 受 控 对 象 为 一 阶 惯 性 过 程 (K2 = 0), 该 控制 器 无 
微分 控制 项 . 
包括 控制 器 (2) 式 , 滤波 器 (3) 式 和 PID 控制 参数 计算 (13) 式 的 这 组 关系 即 为 一 、 二 阶 惯性 加 时 滞 环 节 受 控 对 
象 的 二 自由 度 PID 内 模 控 制 〈IMC-2DOF-PID ) 系统 实现 方案 . 依 (4) 式 得 该 控制 系统 对 激励 输入 和 干扰 的 全 
响应 为 : 


(13) 


5S(X2s+2X2+7) 
Y¥ (5) 和 S (M3S+2N2+7)+e-"s (9) (14) 


(和 2S 十 1) 。 e 一 73 
十 (Na5H ET 


分 析 知 该 响应 除 时 滞 时 间 r 外 与 对 象 的 其 它 参 数 无 关 , 这 就 可 以 引入 一 种 新 的 优化 机 制 : 不 在 现场 针对 具体 
对 象 , 而 是 离线 地 对 一 个 特定 系统 仿真 并 优化 其 响应 ; 通过 制 成 优化 响应 的 模板 , 在 现场 快速 地 辨识 得 实际 对 象 
参数 , 即 可 实现 IMC-2DOF-PID 系统 的 优化 . 原理 见 下 节 . 

3 引入 模范 系统 概念 并 建立 实现 优化 控制 的 新 机 制 
依据 设计 的 IMC-2DOF-PID 系 统 响应 (14) 式 知 , 其 仅 与 滤波 参数 和 实际 对 象 的 时 滞 时 间 7 相关 . 如 令 : 


和 1 LT 入 2 一 Q27 (15) 
有 旦 只 考虑 Y(5) 中 输入 信号 有 R(5) 的 激励 部 分 六 (5S) 得 : 
_ (a2TS+1) e-'SR(S) 
= (aa7S+I)  S(a372S+2027+7)+e-rs (16) 


为 消除 其 中 仍 与 实际 对 象 相关 的 参数 ， 将 时 滞 时 间 7 归 一 化 , 即 得 新 的 系统 响应 Yo(9): 


_ (a25+1). e-SR(S) 
Yor(9) (a1S+1) S(a3S+2a2t+1)te-s (17) 


这 是 个 与 受 控 对 象 参 数 完全 无 关 的 特定 系统 ， 命 名 为 模范 系统 (Modle System). 再 定义 模范 对 象 G。(S) 和 模范 
滤波 器 Gy,(5) 为 : 


—S 


GH) 一 mia (18) 
Grol 0) = et (19) 


就 有 模范 系统 的 结构 如 图 4: 


图 4. 模范 系统 结构 图 


Fig. 4. Model System structure diagram 


当 7( 力 为 单位 阶 距 信号“ 即 R(5) = 93 一) 时 ,(17) 式 的 模范 系统 响应 加-(5) 与 (16) 式 的 内 模 控 制 系统 响 


应 7.(S) 存在 如 下 关系 : 


Yor(T5)= ¥(S)/7 (20) 
即 有 : 
Y(S) = 7Y(7S) (21) 
根据 拉 普 拉 斯 变换 的 相似 定理 有 : 
yr(t) = Yor(t/7) (22) 


该 式 表 明 : 用 实际 对 象 的 时 滞 时 间 对 模范 系统 响应 曲线 做 相似 变换 就 能 得 到 前 述 IMC-2DOF-PID 系统 的 响应 曲 
线 , 既 两 个 系统 的 优化 响应 也 是 相似 的 . 由 此 可 事先 离线 地 仿真 与 具体 对 象 无 关 的 模范 系统 响应 , 按 所 选 控制 性 
能 指标 寻 优 得 到 模范 系统 的 优化 参数 oj , oa. 而 在 获得 对 象 实际 参数 后 , 将 这 组 优化 参数 乘 以 具体 对 象 的 时 滞 时 
间 参 数 7, 可 得 实际 的 内 模 控 制 系统 优化 参数 : 
和 一 adT, NM 一 oT (23) 

将 参数 入 ,入 与 对 象 实 际 参数 Ks, Ki, Ko 一 并 代入 式 (13), 式 (2) 与 式 (3), 所 得 控制 器 和 滤波 器 即 可 构成 优化 的 
内 模 控 制 系统 〈 结 构 见 图 1) . 该 系统 的 优化 响应 为 : 

入 (9) 一 TYorr9) 

时 CS Re es —_R(S) (24) 

(XS+1) S(M25+2X8+7)+e-"rS 

这 就 借助 模范 系统 建立 了 一 种 IMC-2DOF-PID 系 统 特别 的 优化 机 制 : 先 优 化 模范 系统 、 后 辨识 实际 对 象 参 
数 ; 在 完成 对 象 辨 识 时 也 就 实现 了 实际 控制 系统 的 优化 . 

分 析 可 知 模范 系统 的 优化 参数 有 两 个 特点 , 其 一 : 其 不 受 具 体 对 象限 制 而 具有 通用 性 . 其 二 : 允许 有 多 对 不 
同 的 优化 参数 ， 各 自 对 应 一 种 控制 性 能 . 换 句 话说 就 是 第 一 : 优化 时 不 必 考 虑 以 后 控制 的 是 什么 对 象 〈 限 
于 FOPDT、SOPDT 受 控 对 象 范围 ) ; 第 二 ， 辩 识 得 到 的 对 象 参数 配 上 不 同 的 优化 参数 可 直接 实现 不 同 的 控制 
特性 . 这 是 引入 该 特定 系统 的 目的 , 也 是 称 其 为 模范 系统 的 理由 . 默认 模范 系统 的 参数 al, as 是 已 按 某 种 控制 
性 能 指标 优化 的 一 对 值 , 也 称 模范 参数 〈 为 之 后 推导 方便 可 免 去 * 标 记 ) . 

在 举例 说 明 模 范 系 统 的 优化 过 程 前 ， 先 列 出 图 4 模范 系统 的 完整 响应 式 : 

Y,(8) = Yo (3) + Yes) 
一 3 . slogsracati) ros le) 
S(a2S+202+1) (25) 
S(a3S+2a2t+1)te-s (9) 
= Ty(S)R(S) + $1(S)D(S) 
式 中 (35) 为 设 定 值 跟随 特性 ; Ya(5) 为 干扰 抑制 特性 ; 而 Sy(5) 为 灵敏 度 函 数 ; Ty(5) 为 互补 灵敏 度 函 数 . 

一 种 优化 模范 系统 的 方法 是 通过 MATLAB 软 件 按 图 5 搭建 模范 系统 和 计算 TTAE 指标 的 仿真 结构 . 优化 时 
可 先 调整 ws 优化 干扰 抑制 性 能 , 再 调整 al 优化 跟随 性 能 . 然后 通过 此 时 模范 系统 的 开 环 传递 函数 Nyquis 曲线 
到 临界 点 (一 1, 70) 最 短 距离 的 倒数 获得 最 大 灵敏 度 (jMs); 参照 文献 [16]P23 的 算式 求 得 相应 的 幅 值 裕 度 (hg) 与 
相位 裕 度 (7), 评估 系统 的 鲁 棒 性 能 . 

仿真 第 一 步 : 激励 信号 为 零 , 加 阶 跃 干扰 信号 , 调 参数 as 使 受 干扰 影响 的 输出 响应 TTAE 值 最 小 ; 第 二 步 : 
关 断 干扰 源 , 加 阶 跃 激励 信号 , 调整 参数 ai 使 输出 响应 与 激励 信号 之 差 的 [TAE 值 最 小 . 
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图 5. 模范 系统 ITAE 指标 寻 优 仿真 结构 图 


Fig. 5. Model System IIAE index optimization simulation Structure diagram 


寻 优 后 得 到 一 对 模范 参数 及 对 应 的 鲁 棒 性 能 指标 : 
aa = 0.00596, as = 0.383; Jrrag|,,, = 2.871397 
(a2 = 0.1467,2as + 1 = 1.766) (26) 
M, = 1.795; h, > 2.258,y > 32.348° 
优化 性 能 可 有 多 种 选择 , 比如 还 可 改变 仿真 结构 按 1STAE 指标 寻 优 得 到 第 二 对 模范 参数 及 对 应 鲁 棒 性 能 指标 : 
al = 0.00597, as = 0.488; Jrsrag|,,, = 8.057236 
(o3 = 0.2381,2a' +1= 1.976) (27) 
M, = 1.700; h, > 2.429,7 > 34.209° 
或 者 按 77 4 五 寻 优 方式 调整 参数 aa, 按 超 调 指标 为 0% 调整 参数 ai 得 到 第 三 对 模范 参数 及 对 应 鲁 棒 性 能 指标 : 
al = 1.505, az = 0.383; J,_o%|,;, = 4.833082 
(az = 0.1467, 20 + 1 = 1.766) (28) 
M, = 1.795; h, > 2.258,7”> 32.348° 
以 上 三 对 模范 参数 构成 的 模范 系统 均 具 有 一 定 的 鲁 棒 稳 定性 ， 在 后 面 实例 中 要 分 别 用 这 些 模范 参数 实现 对 象 
辨识 与 系统 优化 . 
新 方法 并 不 限制 模范 系统 采用 何 种 手段 优化 , 而 是 根据 实际 控制 需要 选取 . 如 为 保证 鲁 棒 稳定 性 还 可 根据 给 
定 的 系统 最 大 灵敏 度 值 Ms 求 模范 参数 az, 再 优化 模范 参数 ai. 甚至 还 可 用 智能 算法 对 模范 系统 的 参数 进行 优 
化 . 
优化 模范 系统 得 到 模范 参数 只 是 新 方法 的 前 期 工作 , 要 真正 实现 该 IMC-2DOF-PID 系统 的 优化 控制 必须 获 
得 对 象 实际 参数 . 辨识 对 象 实际 参 数 的 方法 将 由 下 节 导 
4 推导 基于 模范 系统 的 对 象 参 数 辨识 算法 
辨识 前 假设 对 象 参 数 初 值 为 Ks, Ki, Ko,, 则 由 模范 参数 可 得 优化 参数 初 值 : Ai(= aiF), 和 (= ao25). 这 组 
初 值 代入 (13), (3), (2) 式 后 得 到 滤波 器 与 控制 器 后 , 图 1 系统 的 信号 激励 响应 为 了 (5): 


7 _ M2S+1 

F(S) = PH x 
K2S?+KIS+Ro i 
X252+(2X2+7)S K2S?+K1S+Ko R(S (29) 
十 Kz2s +K1S+Ko er 局 


对 5?+(2Xo+F)S R25S3+R1ISTRo 
所 选 初 值 要 满足 该 响应 结果 稳定 且 无 偏 ( 即 激励 信号 升幅 后, 经 一 段 时 间 了 有 y(T) = 区 co) = 万 ). 
再 由 模范 系统 获得 响应 模板 Yi (5S) 如 下 : 
(8) = HH. grr jr BS) (30) 
消除 (29), (30) 式 中 激励 信号 有 R(S) 后 合并 得 这 两 种 响应 的 关系 如 下 ( 恒 等 关 系 表明 同 激励 信号 下 必 有 此 式 ): 
| [3352+(2X2+7)S](K2S?+K1S+Ko) | elr-7)sy (8) 


+(K2S?+R1S+Ro)e "S 


一 FS 


(31) 


[X352+(2X2+7)Ste "|(K2S?+K1S+Ko) Yi(S5) 
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其 中 初 值 K2, Ki1, Ko,7， 和 已 知 ;响应 了 (5) 今 
Ki, Ko,T 的 辨识 算法 过 程 如 下 : 

先 令 入 = 了 一 为 时 洁 时 间 偏差 其 可 通过 
重 按 9 各 次 至 归 类 调整 (31) 式 后 得 : 


K2S?|Yor(S)-—Y(S)]e (+)S+ 


K1s|Yor(S)-Y(S)]e (+A)s+ 
Ko[Yor(S)-—Y(S)|]e +)S 
(K2 一 Ka)ss7(s)-R2ss [Yor(S)e-Ss-Y(S)]+ 
(KG 一 Ki)s??(s)-Ris? [Yor(s)e-Ss-Y(S)|+ 
(Ko 一 Ko)sy(s)-Ros|Yor(S)e-Ss—Y(s))] 


步骤 1 式 (32) 中 5S 取 极 限 零 可 消 S253.94 项 得 : 
(Ko = Ko) SY(S)-KoS[Yor(S)e 人 As- 立 (5)|= 
lim Ki+ 区 
S70 2X Eo | [s¥or(s)-s¥(s)] €-(7+A)S 


Ss 
日 拉 普 拉 斯 变换 的 终 值 定理 和 拉 普 拉 斯 反 变 换 得 : 


(Ko — Ko)y ,sn =Ro ly (tA)-9(D) sr 


本 


t—+—A)—y(t—7— 人 A)] 


卫生 [go 人 
一 
22 | yor (t—7—A) -9-7—A)ldt 


tt 之 下 


式 中 了 是 系统 稳定 的 时 间 , 此 时 各 响应 为 RR, 所 以 等 式 右 侧 第 一 、 二 项 为 零 , 第 三 


的 积分 值 不 变 〈 消 除 的 目的 是 便于 实时 计算 ， 


dKo = (Ko = Ko) = RN { [gor (©) g(t)]at 
步骤 2 式 (32) 乘 5-! 并 取 lim 极限 后 , 同 理 得 : 
dK = 一 守 一 Ry (bdt 十 
Ko El (t— 人 A 人) T(t)]dtt Rr) * 


[yor 
是， 
Ki f [yor (lt 下 A) y(t T 人 A)]dt+ 
TT 
Ko ff lyor(t—7—A) ot-7 Ade 


消除 式 右 侧 积分 项 中 偏 移 时 间 , 结合 (35) 式 整理 后 得 


入 
7 dK0o 1 
adKi AKo- 千 /7 (Da RTX 
二 X2 工 
Ko f [yor(t)—y(t)]at+ 
0 
T 


1 

只 

oc, 
KY 
o 
3 
全 
于 
二 


-0a Ro f /I yor (t)—y(t)]dt? 


其 中 dK6o 可 由 (35) 式 计算 . 
步骤 3 式 (32) 乘 9-2 并 取 lim 极限 后 , 同 理 得 : 


-2 -2 更 
dK, 一 ->221 (证 Ha ) f y(t)dtt 
0 


2 和 2 十 元 2 和 2 十 元 


[yor (tA) gat ago f roat2+ 
0 0 


wl 
1 品 全 昌 


[lw 
ba 
[= 
I 


五 
直 Ee +R:) { [Yor (t—A)—y(t)]at+ CAT x 
IT 
K2 [ [yor(t—7—A)—y(t—7—A)]dt+ 
Tr 
Ki [ f [yor(t—7—A)—y(t—7—A)]dt?+ 
TT 
Ko 厂矿 厂区 (一 = A)— y(t— 7—A)]ats 
000 


实测 系统 响应 曲线 与 响应 模板 曲线 后 求 得 . 为 后 续 逐 


三 (X35+2X2+7) x 


步 推 


项 中 消除 偏 移 时 间 (F 1 


g(t) 可 实测 :模板 二 (5) 今 yo (t) 可 解 算 ; 对 象 实 际 参数 K， 


导 方 便 ， 


(32) 


(33) 


下 同 ), 整理 (34) 式 后 即 得 对 象 参数 Ko 与 初 值 Ko 的 差 值 


(36) 


(37) 


(38) 
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消除 式 右 侧 积分 项 中 偏 移 时 间 , 结合 (35), (37) 式 整理 后 得 : 


x2 
i yo A2 法 
aK, (an Ko+ARi+ ao) FR X 
2X2 


[at (232—A) FitRBRot © Ro] f lor() -0 


+[Ki A)Eol 了 1 gor (t)—g(t)]dt? (39) 


+(2X2— 
e 

+Ko ff ff [yor(t)—y(t)]ats 
5 0 


~ x2aKk 更 TT 
一 去 |AKo Out+( Re Ha ) f yt)adttdaKo ff y(t)at 
0 0 0 0 


其 中 dKo, dK 可 由 (35), (37) 式 计算 . 
至 此 可 在 同一 信号 aa dale mal 曲线 , 同时 测算 入 、 系 统 响 应 与 
响应 模板 的 积分 和 重 积分 ; 再 经 (35), (37), (39) 式 解 算出 实际 对 象 参数 的 各 差 值 , 并 由 下 式 算 出 对 象 各 参数 : 


7T 一 示 十 人 

Ko = Ko+ dkKo (40) 
Ki= Ki+daKi 

K, = K, + dkK, 


从 理论 上 说 一 次 计算 就 可 得 到 实际 对 象 参 数 , 但 考虑 实际 系统 存在 测量 噪声 或 干扰 , 解 算 的 对 象 参数 值 也 有 
误差 , 故 将 其 作为 新 一 轮 初 值 迭 代入 循环 的 辨识 过 程 , 并 引入 模范 响应 相似 偏差 评价 值 (简称 相似 差 值 ) J (TT)， 
见 (41) 算式 , 以 其 值 最 小 作为 循环 结束 的 条 件 : 

J(T) =/ [907 (t) 9()dt (41) 
循环 结束 后 将 所 获得 的 对 象 参数 r, Ko, Ki, Kz 代入 (15), (13) 式 , 即 可 获得 优化 的 滤波 器 和 控制 器 参数 , 同时 使 
空 制 系统 实现 优化 响应 . 这 就 构建 了 一 种 基于 模范 系统 同步 实现 受 控 对 象 参数 辨识 及 其 控制 系统 优化 的 新 方法 . 


该 方法 在 现场 的 迭代 辨识 过 程 与 已 有 优化 方法 中 的 极 值 搜索 (Extremum-Seeking) 过 程 看 似 很 像 , 但 实际 上 
有 本 质 不 同 : 前 者 参照 的 是 已 优化 的 响应 模板 , 而 后 者 是 针对 具体 对 象 控制 系统 的 某 个 性 能 指标 . 
5 实例 仿真 及 分 析 


为 验证 展示 上 述 机 理 及 算法 的 可 行 性 , 本 节 借 助 MATLAB 软件 仿真 某 对 象 、 搭 建 IMC- 2DOF- PID 控制 系 
统 及 模范 系统 相似 的 响应 模板 , 结构 见 图 6; 建立 图 7 辨识 算法 模块 , 该 模块 左 侧 输 入 模范 参数 和 对 象 初 值 , 左下 
输入 控制 系统 和 响应 模板 的 输出 , 右 侧 输 出 对 象 参数 差 值 和 相似 差 值 . 通过 仿真 信号 激励 系统 和 循环 数据 迭代 
过 程 , 最 终 辨识 出 对 象 参 数 并 实现 系统 的 优化 控制 . 详 见 以 下 实例 : 

辨识 时 响应 模板 采用 (26) 式 的 模范 参数 ; 控制 系统 中 加 入 约 2% 幅 值 白 噪声 模拟 扰动 。 


Dis-y(t )&yor(t) 


图 6. 控制 系统 和 响应 模板 结构 


Fig. 6. Control system and the response template structure 
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图 7. 辨识 算法 模块 


Fig. 7. Identification algorithm module 


例 1 某 二 阶 惯性 时 延 对 象 为 : 
Gi(O) = ma (42) 
对 象 实际 参数 : 7 = 9.0; Ko = 6.0; Ki = 13.0; Kz ==7.0; 无 死 区 . 辨识 前 设 对 象 参 数 初 值 : 7 = 6.0; 
Ko = 3.0; Ki = 12.0; Ks = 0.0; 构成 初 态 过 阻尼 系统 . 经 多 轮 仿真 迭代 后 的 辨识 数据 见 表 1. 每 轮 对 象 参数 初 
直 经 差 值 修正 后 作为 下 轮 初 值 重 新 辨识 , 直至 确认 得 到 最 小 相似 差 值 后 结束 . 


表 1 初 态 过 阻尼 系统 辨识 数据 
Table 1 Identification data of initial over-damped system 


~ 


对 象 参数 T K0 Kl K2 | 相似 差 值 J 
(实际 值 ) (9.00) (6.0000) (13.0000) (7.0000) | (1.523849) 
首 轮 | 初 值 6.00 3.0000 12.0000 0.0000 | 11.737854 
差 值 | 3.051 2.99992 0.74194 4.19535 V 
二 轮 | 初 值 9.05 5.9999 12.7419 4.1954 | 1.538553 
差 值 | 0.073 0.00009 -0.45513 1.18036 V 
三 轮 | 初 值 9.12 6.0000 12.2868 5.3758 | 1.523964 
差 值 | -0.060 0.00000 ”0.38332 0.85334 V 
四 轮 | 初 值 9.06 6.0000 12.6701 6.2291 | 1.523772 
差 值 | 0.011 0.00000 -0.07076 -0.13267 V 
五 轮 | 初 值 9.07 6.0000 12.5993 6.0964 | 1.523791 
差 值 - - . 区 
辨识 结果 9.06 6.0000 12.6701 6.2291 | 1.523772 
表 中 符号 


V 表 示 此 轮 的 相似 差 值 仍 最 小 ,继续 辨识 ; 
X 表 示 此 轮 的 相似 差 值 非 最 小 , 结束 辨识 ; 
- 表示 参数 修正 差 值 无 效 .( 下 表 同 ) 


202211.00245Vv2 


chinaXiv 


经 五 轮 辩 识 过 程 , 获得 对 象 参数 : 了 = 9.06; Ko = 6.0000; Ki 


从 开始 的 11.737854 减 到 最 小 值 1.523772 ( 受 噪 声 影响 不 为 零 , 仿真 可 知 对 象 实际 参数 的 ,jy 值 1.523849, 见 表 
内 J 栏 括号 中 值 ). 优化 前 后 曲线 变化 反映 在 图 8 中 , 黑 线 为 模板 响应 , 灰 线 为 控制 系统 
8(b) 表明 第 二 轮 经 首次 优化 后 控制 系统 响应 与 模板 响应 基本 重合 , 且 与 图 8(c) 最 后 


优化 的 响应 曲线 相 比 难以 看 出 明显 差异 , 只 是 后 者 J 值 (1.523772) 小 于 前 者 的 值 (1.538553). 


值 使 系统 初 态 为 过 阻尼 ; 图 


经 七 轮 辨识 过 程 ， 获 得 对 象 参 数 : + = 9.04; Ko = 6.0000; Ki = 12.7965; Ks = 6.5201; 相似 差 值 J 


二 12.6701; K。 = 6.2291. 各 轮 相似 差 值 J 


(a) 优化 前 的 响应 


(b) 首次 优化 后 的 响应 


: : l . B 
0 50 100 150 200 250 300 


(0) 最 终 优化 的 响应 
图 8. 初 态 过 阻尼 系统 辨识 优化 


响应 


Fig. 8. Response in identification optimization of initial over-damped System 


如 选 参 数 初 值 : + = 6.0; Ko = 8.0; Kj = 19.0; Ks = 0.0; 则 系统 初 态 为 欠 阻 尼 , 也 能 完成 辨识 和 优化 . 几 
轮 的 辨识 数据 见 表 2. 


向 应 . 图 8(a) 显示 首 轮 初 


从 29.217253 减 至 最 小 值 1.523792. 初 态 欠 阻尼 系统 优化 前 后 变化 如 图 9 所 示 , 图 9(b) 中 首次 优化 后 系统 响应 就 


基本 与 响应 模板 
一 致 . 说 明 图 8(c)， 


图 9(c) 都 实现 了 优化 (模范 ) 的 响应 . 


合 . 虽然 初 态 两 种 阻尼 使 系统 响应 差异 极 大 , 参照 模板 辨识 使 表 1. 表 2. 的 两 组 相似 差 值 趋 


表 2 初 态 欠 阻尼 系统 辨识 数据 


10 


近 
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Table 2 Identification data of initial under-damped system 


对 象 参数 T K0 KI K2 | 相似 差 值 J 
(实际 值 ) (9.00) (6.0000) (13.0000) (7.0000) | (1.523849) 
首 轮 | 初 值 6.00 8.0000 19.0000 0.0000 | 29.217253 
差 值 | 2.949 -2.00536 -4.85624 -0.85487 V 
二 轮 | 初 值 8.95 5.9946 14.1438 0.0000 | 1.766383 
差 值 | 0.257 0.00539 -2.36195 4.78363 V 
三 轮 | 初 值 9.21 6.0000 11.7818 4.7836 | 1.524232 
差 值 | -0.065 0.00000 ”0.38033 0.43931 V 
四 轮 | 初 值 9.14 6.0000 12.1621 5.2229 | 1.524015 
差 值 | -0.091 0.00000 0.57873 1.15835 V 
五 轮 | 初 值 9.05 6.0000 12.7408 6.3813 | 1.523806 
差 值 | -0.011 0.00000 0.05565 0.13884 V 
六 轮 | 初 值 9.04 6.0000 12.7965 6.5201 | 1.523792 
差 值 | 0.002 0.00000 -0.02361 -0.05205 V 
七 轮 | 初 值 9.04 6.0000 12.7729 6.4680 | 1.523796 
差 值 - - - 3 X 
辨识 结果 9.04 6.0000 12.7965 6.5201 | 1.523792 


05 人 hn 看 1 | 
0 50 100 150 200 250 300 
(a) 优化 前 的 响应 


(b) 首次 优化 后 的 响应 


Se 


(c) 最 终 优 化 的 响应 
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图 9. 初 态 欠 阻尼 系统 辨识 优化 中 响应 


Fig. 9. Response in identification optimization of initial under-damped system 


在 已 获 对 象 参 数 的 基础 上 , 模范 参数 换 成 (27) 式 中 ISTAE 优化 指标 的 一 对 , 可 立即 得 到 图 10(a) 中 新 优化 

案 的 系统 响应 . 因 优 化 指标 不 同 , 响应 动态 指标 明显 改变 . ITAE 的 超 调 量 14.3%, 上 升 时 间 Tr = 24s; 

0 的 超 调 量 10.5%, 上 升 时 间 Tr = 26s. 如 换 成 (28) 式 中 一 对 模范 参数 , 则 可 得 零 超 调 量 优化 方案 的 系 
统 响应 , 见 图 10(b) 中 曲线 . 其 超 调 量 0%, 上 升 时 间 Tr = 54s. 


0 50 100 150 200 250 300 


(a) 按 1STAE 指标 优化 的 系统 响应 


0 50 100 150 200 250 300 


(b) 按 超 调 量 为 0% 优化 的 系统 响应 
图 10. 线性 对 象 不 同 优化 方案 的 系统 响应 


Fig. 10. System response of linear object with different optimization schemes 


这 说 明 借助 模范 系统 的 模范 参数 可 快速 改变 IMC-2DOF-PID 系统 的 优化 特性 . 
例 2 局 部 含 非 线性 的 高 阶 对 象 为 : 


a 0.5S+1 一 9.05 
Gpaz(9) = 555T7.055 二 13.05T6.05 (43) 


对 象 实际 参数 :7 =9.0; Ko = 6.0; Ki = 13.0; Kz==7.0; 加 上 零点 (0.55 十 1), 高 阶 项 (1.553) 和 死 
土 0.6s). 设 初 值 : + = 6.0; Ko = 3.0; Ki = 12.0; Ks = 0.0, 该 对 象 的 辨识 数据 见 表 3. 

经 七 轮 辨 识 过 程 , 获得 等 效 的 对 象 参 数 : 了 = 9.78; Ko = 6.6000; Ki = 10.5559; KK。 = 36.2803. 优化 前 后 
响应 曲线 变化 如 图 11 所 示 . 受 非 线性 等 因素 影响 首次 优化 的 响应 与 响应 模板 有 部 分 明显 偏差 , 但 继续 至 最 终 优 
化 后 也 实现 了 重合 , 只 是 相似 差 值 稍 大 些 .通过 变换 模范 参数 仍然 可 实现 不 同 特性 的 系统 响应 , 见 图 12. 


区 


> 


表 3 局 部 含 非 线 性 高 阶 对 象 辨识 数据 
Table 3 Nonlinear high-order object identification data 
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对 象 参数 T K0 K1 K2 | 相似 差 值 J 
(实际 值 ) (9.00) (6.0000) 〈13.0000) (7.0000) (-) 
首 轮 | 初 值 6.00 -3.0000 12.0000 ”0.0000 | 13.822837 
差 值 | 4.401 3.59990 -2.25181 -0.28465 V 
二 轮 | 初 值 10.40 6.5999 9.7482  *0.0000 | 2.660447 
差 值 | 0.499 0.00010 -7.48564 18.76837 V 
三 轮 | 初 值 10.90 6.6000 2.2626 18.7684 | 2.059064 
差 值 | 0.237 0.00000 -0.88889 17.51192 V 
四 轮 | 初 值 11.14 € 6.6000 1.3737 36.2803 | 1.919147 
差 值 | -0.534 0.00000 4.11257 11.36436 V 
五 轮 | 初 值 10.61 6.6000 5.4863 *36.2803 | 1.736454 
差 值 | -0.826 0.00000 ”5.06949 8.73029 V 
六 轮 | 初 值 9.78 6.6000 10.5559 *36.2803 | 1.694844 
差 值 | -0.307 0.00000 1.43994 2.88712 V 
七 轮 | 初 值 9.47 6.6000 11.9958 39.1674 | 1.752308 
差 值 - - = = X 
辨识 结果 9.78 6.6000 10.5559 36.2803 | 1.694844 
* 必 要 时 可 对 参数 限定 取 值 范围 , 如 : 0 < K2 < 40 


(0) 最 终 优化 的 响应 
图 11. 局 部 含 非 线性 对 象 系统 辨识 优化 中 响应 


Fig. 11. Response in System identification Optimization with nonlinear high-order object 
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(a) 按 757 4 已 指标 优化 的 系统 响应 


50 100 


150 200 250 300 


(b) 按 超 调 量 为 0% 优化 的 系统 响应 


图 12. 非 线性 高 阶 对 象 不 同 优化 方案 的 系统 响应 


Fig. 12. System response of nonlinear high-order object with different optimization schemes 


从 以 上 变换 后 的 控制 效果 看 , 即使 是 局 部 含 非 线性 、 零 点 及 高 阶 环节 的 对 象 也 可 等 效 为 二 阶 惯性 时 消 环 节 


模型 . 
以 上 几 种 情况 只 经 过 


、 两 次 优化 后 , 系统 响应 


昌 线 与 优化 响应 模板 就 基本 重合 ; 总 次 数 也 都 在 个 位 数 以 


内 , 相 比 极 值 搜索 动 邱 几 十 其 至 上 百 的 循环 次 数 , 明显 快速 高 效 . 说 明 前 期 对 模范 系统 的 优化 过 程 , 已 将 系统 优 
化 工作 与 现场 对 象 辨识 工作 脱离 ; 虽然 看 似 辨识 与 优化 最 后 同时 完成 , 实际 是 将 前 期 优化 模范 系统 得 到 的 响应 


模板 作为 系统 响应 的 “模具 


同时 也 实现 了 优化 的 内 模 
6 总 结 


99 


“范本 ”进行 辨识 


空 制 ; 工作 效率 因此 提高 . 


与 整定 , 当 系统 响应 和 模板 趋 于 一 致 时 得 到 对 象 实际 参数 ， 


新 方法 不 同 于 PID 控 制 或 内 模 控 于 


| 传统 的 实现 流程 与 整定 方式 , 而 是 用 模范 系统 这 个 特定 模型 离线 做 优化 时 


得 到 的 优化 响应 模板 和 相关 算法 引导 ， 现 场 在 线 快 速 地 辨识 对 象 并 同时 整定 系统 优化 控制 参数 . 优点 是 : 1. 模 


范 系统 优化 方法 不 受 实际 受 控 对 象限 制 ; 2. 离线 优化 模范 系统 方式 易于 实施 又 不 会 影响 生产 安全 ; 3. 更 换 其 它 
受 控 对 象 不 必 重 新 优化 模范 系统 ; 4. 在 响应 模板 和 算法 引导 下 在 线 辨 识 与 整定 易于 实施 且 快 捷 高 效 ， 也 便于 控 


制 系统 进行 后 期 维护 ，5. 辨识 对 象 参 数 后 换 不 同 优化 模板 可 直接 实现 新 控制 性 能 (不 必 重 新 辨识 整定 ). 此 法 中 


的 模范 系统 适用 于 工业 生产 过 程 中 多 数 常 见 的 二 阶 惯性 加 时 灌 环 节 或 局 部 含 非 线 性 、 零 点 高 阶 环节 受 控 对 和 象 . 


该 方法 无 需 精确 对 象 参数 , 参数 整定 以 算法 为 依据 , 降低 了 对 现场 实施 人 员 的 技术 与 经 验 要 求 , 减少 人 为 误 操 作 
影响 . 其 控制 器 的 内 模 控制 特性 ， 


二 


统 软 、 硬 件 实现 , 便于 在 生产 中 


象 ， 该 方法 应 采用 预 补偿 处 至 
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